Trace anomaly, thermal power corrections 
and dimension two condensates in the deconfined phase 



o 
o 



E. Megi'afB 

Nuclear Theory Group, Physics Department, Brookhaven National Laboratory, Upton, New York 11973 USA 

E. Ruiz Arriolsl!! and L.L. SalcedcS 

Departamento de Fisica Atomica, Molecular y Nuclear, 
Universidad de Granada, E-18071 Granada, Spain 
(Dated: September 10, 2009) 

The trace anomaly of gluodynamics on the lattice shows clear fingerprints of a dimension two 
condensate above the phase transition. The condensate manifests itself through even powers of the 
inverse temperature while the total perturbative contribution corresponds to a mild temperature 
dependence and turns out to be compatible with zero within errors. We try several resummation 
methods based on a renormalization group improvement. The trace anomaly data are analyzed and 
compared with other determinations of the dimension two condensate based on the Polyakov loop 
and the heavy qq free energy, yielding roughly similar numerical values. The role of glueballs near 
the transition is also discussed. 
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I. INTRODUCTION 



The physics of the Quark-Gluon Plasma has turned 
out to be more difficult than initiaUy expected (for a re- 
cent review see e.g. [l| and references therein). While 
it has been noted 0] that the failure of the perturba- 
tive description could have been anticipated on general 
grounds, it has not been obvious what is the form and 
furthermore the dynamical origin of such a strongly non- 
perturbative behaviour. On the other hand it is by now 
clear from lattice calculations at finite temperature (for a 
review on recent developments see e.g. [3] and references 
therein) that while in QCD there is a cross-over, in pure 
Yang-Mills theory or gluodynamics for iVc > 3 there is 
a ffi'st order phase transition. Actually, in the simpler 
case of gluodynamics the deconfinement phase transition 
is monitored by the expectation value of the Polyakov 
loop which acts as an order parameter of the associated 
breaking of the center 1^{Nc) symmetry. 

The finding of power corrections above the phase tran- 
sition [4] and the possible explanation in terms of a 
dimension-2 condensate of the dimensionally reduced 
theory of such a special object suggests pursuing the same 
idea in other thermodynamic quantities. Power correc- 
tions are ubiquitous at high energies and zero temper- 
ature through condensates being the remnant of non- 
perturbative effects. From this viewpoint there is no 
reason why they should not be present at high tem- 
peratures as genuinely non-pcrturbative effects. In the 
present paper we detect and analyze power corrections 
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on the trace anomaly density at finite temperature as 
obtained in lattice calculations, and which can only be 
effectively explained by a dimension two condensate. Un- 
like the Polyakov loop, the trace anomaly corresponds to 
a physical and measurable quantity not only on the lat- 
tice but also experimentally in the QCD case. In this 
way we extend our study of the non-perturbative effects 
above the phase transition started for the Polyakov loop 
and the qq free energy 0, d, [B| • A brief account of our 
results have been presented in Refs. 0, Si- 

The paper is organized as follows. In Section |TT] we 
provide mounting evidence on the dominance of power 
corrections for the trace anomaly using available lattice 
data for gluodynamics 0. The thermodynamics of the 
phase transition is analyzed in Section IIIII following an 
interesting insight by Pisarski flO]. The physics below 
the phase transition in the confined region should be de- 
scribable in terms of colour singlet glueball states, an is- 
sue that we address within a large Nc perspective where 
the interactions among glueballs are much suppressed. 
In passing we note a Hagedorn looking glueball exponen- 
tially growing spectrum with practically no impact on 
the confined region. We analyze in Section IIVI our re- 
sults from the point of view of analytical schemes based 
on perturbation theory. Actually, a common deficiency 
of the poorly converging perturbation theory is the oc- 
currence of unpleasant infrared divergences and manifest 
spurious scale dependence. Perturbative resummations 
such as Hard Thermal Loop remove the infrared singular- 
ities due to Debye screening but provide too weak a signal 
as compared to the lattice data. Our analysis suggests 
that those methods might in fact describe the deconfined 
phase after the power corrections have been subtracted 
off. To reinforce this conclusion a renormalization group 
improvement of the perturbative free energy is under- 
taken in this Section. We find that besides implement- 
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ing renormalization scale independence, this resumma- 
tion mimics the Hard Thermal Loop calculations within 
estimated uncertainties, and as a consequence similar 
conclusion follows: power corrections dominate the de- 
confined phase while perturbatively based physics plays 
a minor role wherever it may be considered applicable. 
The dimension two condensate has naturally explained 
the presence of power corrections for the Polyakov loop 
and singlet qq free energy 0, @ . The idea is tried out in 
Section |V] for the trace anomaly where we indeed show 
that there is a satisfactory numerical agreement for the 
dimension two condensate with that of the Polyakov loop 
and the qq free energy. Finally in Section IVTl we give our 
main conclusions. In Appendix [A] we rederive the trace 
anomaly at finite temperature for reference. In Appendix 
[B]the possible Hagedorn pattern of the glueball spectrum 
is analyzed. In Appendix [C] we study an improvement 
of the perturbative expansion of the free energy at finite 
temperature based on renormalization group invariance. 



II. THERMAL POWER CORRECTIONS IN 
THE TRACE ANOMALY 

For zero or infinite quark masses (gluodynamics) the 
classical scale invariance of QCD is broken by quantum 
corrections due to the necessary regularization yielding 
a trace anomaly [ll[. In the context of QCD sum rules 
phenomenology the issue has been discussed in Ref. [l3 
(for a review see e.g. [l^l)- At finite temperature the 
trace of the energy momentum tensor is related not only 
to the trace anomaly but also to the difference e — 3P 
with e the eiiergy density and P the pressure [l3l (see 
also Ref. [isl Il6| ) providing a measure of the interaction 
as well as the anomalous breaking of scale invariance. An 
early lattice calculation of both energy and pressure was 
first undertaken one decade ago [§, [l3] (see also more re- 
cent calculations [IM HO) IIH ) exhibiting the expected 
deconfinement phase transition although far still from 
a free gluonic plasma. Finite temperature Ward identi- 
ties and low energy theorems based on the trace anomaly 
have been deduced in Refs. [2^ ,23il in the continuum and 
extended to Euclidean lattices [2J|. Renormalization is- 
sues and the identity of thermodynamic and hydrostatic 
pressure was discussed in Ref. [25| within the real time 
formulation. Phenomenological implications in the con- 
text of resonance gas model are analyzed in [2^ . A good 
knowledge of the trace anomaly is crucial to understand 
the deconfinement process, where the non perturbative 
(NP) nature of low energy QCD seems to play a promi- 
nent role 0. 

The derivation of the well known relation 

e-3P=^((G«j2), (2.1) 

between the trace anomaly and the beta function of glu- 
odynamics is reviewed in Appendix [XI 



The trace anomaly density has been computed at fi- 
nite temperature on the lattice for pure gluodynamics 0] . 
The standard plot of 

e - 3P d f P\ , , 

A^^ = r-(^j (2.2) 

as a function of T is shown in Fig. [TJ Below the critical 
temperature the trace anomaly A is very small though 
not exactly zero. A increases suddenly near and above Tc 
by latent heat of deconfinement, and raises a maximum 
at T K, 1.1 Tc. Then it has a gradual decrease reaching 
zero in the high temperature limit. The high value of A 
for Tc < T < iTc corresponds to a strongly interacting 
Quark-Gluon Plasmapicture. 

In previous works [J, |6| we have detected the presence 
of inverse power corrections on other thermal observables. 
Guided by our previous experience we represent in Fig. [5] 
the dimensionless interaction measure (e — 3P) /T'' as a 
fmiction of (in units of Tc). This plot exposes an 

unmistakable straight line behaviour corresponding to a 
power correction of the form 

^=«a + &a(^)' (2.3) 

in the deconfined phase region slightly above the critical 
temperature. A direct fit of the lattice data {N^ x Nr — 
W X 4) for 1.13 < T/Tc < 5.0 yields 

flA = -0.041(17), 6a = 3.99(5), xVdof = 5.0 (2.4) 

with correlation parameter r(a, 6) = —0.78. Likewise, 
the fit for {N^ x Nr ^ 32^ x 8) and 1.13 < T/Tc < 4.5 
yields 

flA = -0.02(4), 5a = 3.46(13), xVdof = 0.35 (2.5) 

with r(a, b) = —0.73. The fit is more accurate in the sec- 
ond set of lattice data. This set corresponds to a more 
precise determination of the trace anomaly on the lattice. 
In order to estimate the result in the continuum limit, we 
can assume that the difference between the two lattice re- 
sults is entirely due to finite cutoff effects. Assuming fur- 
ther that the corresponding leading effect goes as l/N^, 
yields the estimate qa = -0.02(9) and 5a = 3.28(27) in 
the continuum limit. 

A fit of the data completely excludes the existence of 
the odd power corrections Tc/T and {Tc/Tf in A. We 
have attempted to determine the coefficient of a possi- 
ble quartic correction, appending formula (|2.3p with a 
term ca {Tc/Tf. A fit of lattice data for iV^ = 8 re- 
sults in aA = -0.05(7), 5a = 3.7(7) and ca = -0.3(9) 
with x^/dof — 0.31 and correlations r-(a, 5) — —0.85 , 
r(5, c) = —0.97 and r{a, c) = 0.75. So, more accurate 
data are desirable in order to identify contributions from 
condensates of dimension 4. 

In any case it is clear the existence of power corrections 
which are beyond the scope of the radiative corrections 
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FIG. 1: The trace anomaly density (e — SP)/T'^ as a function 
of T (in units of Tc ). Lattice data are from 9] for A''^ x Nr — 
16^ X 4 and 32^ x 8. The fits using Eq. dO)) are plotted. 2- 
loop result for the trace anomaly density in HTL perturbation 
theory, from [2^, is shown as a shaded band that correspond 
to varying p by a factor of two around 2nT. 



Clearly this expression cannot be used at very high tem- 
peratures unless OA = 0, a result suggested by the fit. On 
the opposite extreme, integrating from very high ener- 
gies, where radiative corrections should dominate, yields 

where 

bo = (iVc - l)27rV45 « 3.51 . (2.8) 

Matching at some temperature T = Tp, gives from 
Pc = P{Tc), bo and Apcrt{T) (assuming ca = 0). Clearly, 
since Pc > 0, and Apcrt(T) > 0, we expect 5a < 6o- 
(Note that the power correction should be negligible at 
the matching temperature.) The result from the fit al- 
most saturates the inequality since the continuum limit 
extrapolation yields 6a = 3.28(27) which suggests that 
both Pc and AportiT) must be very small. The analysis 
of Section Hill confirms these issues. 



III. THERMODYNAMICS OF THE PHASE 
TRANSITION 
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FIG. 2: The trace anomaly density (e — 3P)/T"* as a function 
of (in units of Tc ). Same data as in Fig. [T] 

accounted for in perturbation theory. It should be men- 
tioned that similar power corrections have been identified 
in the trace anomaly in gluodynamics in 2 + 1 dimensions 

Integrating the trace anomaly, as modeled in (|2.3p . 
from T = Tc and using the continuity of the pressure 
across the phase transition yields the following form for 
the pressure 

P(T) bA T2\ , /r\ P{Tc) 

(2.6) 



The analysis of the previous section leaves little doubt 
on the existence of power corrections in gluodynam- 
ics down almost to the critical temperature. Moreover, 
at the maximum temperature measured on the lattice, 
T = 5Tc, the power contribution still dominates the re- 
sult, and both terms become comparable only at higher 
temperatures, T ^ 10 Tc. One sensible question to ask 
is what is the temperature where perturbation theory 
applies. We address this question in more detail in Sec- 
tion Hvl 

If for the moment perturbative corrections are disre- 
garded (they are small anyhow), it is worth reviewing 
why the previous results are completely against the idea 
of a gluon plasma. Actually, the standard textbook ar- 
gument requires the bag pressure to equilibrate the free 
gluon gas at the critical temperature 

P(r) = Pgluons(T) - P, T>Tc, (3.1) 

yielding the trace anomaly 

A-^, T>Tc. (3.2) 

As already shown, this behaviour is excluded by lattice 
data. 

In this section we elaborate on a suggestion by Pis- 
arski (loj . 

P{T) = j{T'^ ~T^T^), (3.3) 

which achieves a fair description of the data assuming 
infinitely heavy glueball masses, negligible radiative per- 
turbative corrections, and a temperature dependent bag 
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energy, called fuzzy bag, 



bo. 



^fuzzy 2 



(3.4) 



As shown in Section |Vl diniension-2 condensates provide 
a natural explanation of this fuzzy bag picture including 
the correct power behaviour and coefficient, while simul- 
taneously account for power corrections in the Polyakov 
loop and the gg-potential at a quantitative level. A study 
of the thermodynamics of this fuzzy constant is also made 
in Ref. |22i] within the holographic QCD approach based 
on the AdS/CFT correspondence. 



A. Infinitely massive glueballs 

At very high temperatures, where radiative correc- 
tions become negligible, the free gluonic gas relation 
PiT) ~ (iV2 - l)7r2T4/45 holds. On the other hand, 
in the confined phase, i.e., below the critical temper- 
ature Tc, the spectrum is believed to be saturated by 
glueballs (for reviews see e.g. Ref. [13, HH). Since the 
pressure must be a continuous function across the phase 
transition, one finds 



P{Tc) — -PglucballsCTc 



(3.5) 



However, the lightest glueball has J^*^ = 0++ and 
Mq++ = 1.73(l)GeV. This is much heavier than the criti- 
cal temperature Tc ~ 270 MeV (a mysteriously disparate 
scale [32|). Thus A is very small (though not exactly 
zero) in the confined phase. Taking the infinite glueball 
mass limit produces Pgiucbaiis = 0. These conditions are 
met by the interpolating functions of the form 



P 



{Nl ~ 1)^^ 
45 



rp4 



T 



T>Tc, (3.6) 



with arbitrary positive n. This parameterization yields 
A = '"'"'^ ( 4^ ) , T>Tc, (3.7) 



45 



T 



which for n = 2 corresponds to take = and 6a — 
3.51 in Eq. (|2.3p . in excellent agreement with the fit to the 
lattice data of the previous section. This thermodynamic 
consistency does not explain however why there is a power 
correction with n — 2. A fit of Eq. (|3.7p to the data 
of Fig. [2] for the same range as in Eq. ((23)) . l.lSTc < 
T < 4.5ATc, yields n = 1.97(5) with a slightly larger 
xVdof = 0.62. 



B. Finite mass glueball spectrum 

Given the goodness of the description we may try to 
improve by actually computing the pressure of an en- 
semble of finite mass glueballs in the confined phase us- 
ing the 12 glueball spectrum currently determined on the 



lattice [33, l3j| . Quite generally such a calculation would 
require the use of the quantum field theoretical version 
of the virial expansion [3^ and thus to take into account 
corrections from binary, ternary and higher order colli- 
sions. For instance, for a binary collision the threshold 
for the first possible lowest glueball-glueball resonance is 
located at twice the glueball mass 2Mq++ ~ 3GeV which 
lies above the 5th glueball and is expected to produce an 
insignificant contribution for T < Tc w 270 MeV. For- 
tunately, the first four low-lying glueballs computed on 
the lattice are kinematically stable bound states. Actu- 
ally, within a large Nc framework one has a gas of sta- 
ble and non-interacting glueballs since Mgb ^ and 
V^b-gb ^ So it seems safe to treat the glueballs as 

a gas of free bosons. This yields 



glueballs 



)3 gEfc/T „ l'' 



(3.8) 



where the index i runs on the glueball species: Mi de- 
notes the glueball mass and gi = 2Ji + 1 the correspond- 
ing angular momentum degeneracy. A convenient low 
temperature expansion gives 



>^luoballs(J ) = 2^ 2^ n 



2n27r2 



V T 



(3.9) 



where the sum over n corresponds to the thermal loops 
and K2 is a modified Bessel function of the second kind. 

Using known spectrum of the lightest 12 glueball 
species 0, [U one can evaluate Pc — -Pgiucbaii(Tc)- This 
requires fixing the critical temperature scale Tc where the 
trace anomaly is evaluated [9| on the one hand and the 
Sommer scale rg defined as TQVqq{rn) = 1.65 [36[, where 
the glueball spectrum is obtained [H, HI], on the other. 
This is necessary since the numerical value of Pc is quite 
sensitive to the glueball mass to critical temperature ra- 
tio, M,/Tc. 

Boyd et al. find TJV^ = 0.629(3) while Luscher 
et al. [33] obtain ro\/a- = 1.22(8). This is consis- 
tent with the elaboration of Teper [30] where r^^/a = 
1.195(10), Tc/^/a = 0.640(15) and Mo++ / ^/^ = 3.52(11) 
are quoted as well as the fixed point action of Niede- 
meyer et al. (ssj simulation where the values VQ^fa = 
1.197(11), Tcl^fa = 0.624(7), tqTc = 0.750(5) and scalar 
roAfo++ = 4.12(21) and tensor roM2++ = 5.96(2^ low- 
lying glueball masses. Along similar lines Necco [33] , gets 
tqTc = 0.7498(50). These latter values agree within er- 
rors with the 12 glueball spectrum simulation [53, 13^ ] 
where roMo++ = 4.21(11) and roA/2++ = 5.85(2) vs 
roMo++ = 4.16(11) and roM2++ = 5.85(5) are obtained 
respectively. Using for definiteness tqTc = 0.75 we get 
the result 



Pc = Pglucball(Tc)= 0.01(1) 



(3.10) 



where the error has been estimated from the uncertainty 
in the glueball masses. This value of the pressure is con- 
sistent with the fit to the data, although the glueball 
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spectrum contribution is not so strongly constrained by 
the trace anomaly data. About 80% of the total contribu- 
tion to OA is given by the lightest scalar O"'"^ and tensor 
2++ glueball states, confirming the marginal role of the 
excited glueball spectrum as well as possible two-glueball 
resonances. 

The trace anomaly due to the glueballs reads 

-i-^-^t^.^x^m- (3.") 

At the lowest lattice temperature T — 0.89 Tc one has 
A < 0.08. This sets an upper bound since all contribu- 
tions in Eq. (|3.1ip are positive. Taking the known glue- 
ball spectrum one gets Agb(0.89rc) = 0.03(6). However, 
the steep raise below Tc is not reproduced by the glue- 
ball spectrum, basically due to the heavy values of the 
masses. In this context it has been suggested [1^ that 
above the phase transition glueballs lower their masses 
providing an enhancement already below the phase tran- 
sition. The behaviour and effect of glueballs in the de- 
confined phase has also been considered in [4l[. Note 
also that the glueball gas formula neglects glueball inter- 
actions which are nominally Assuming that the 
multiplicity in the glueball spectrum is 0{N^) this yields 
a trace anomaly A ~ 0{N^) for T <Tc, while we clearly 
have A = 0{N^) for T > Tc. 

Further considerations on the glueball spectrum and its 
possible Hagedorn structure are presented in Appendix 
IbI In all, our analysis above explains why the pressure 
of finite mass glueballs is compatible with zero within 
errors both from the trace anomaly lattice calculation 
as well as from the direct glueball spectrum estimate, 
and hence supports the infinite mass glueball assumption 
underlying the model of Ref . [13] ■ 

IV. RENORMALIZATION GROUP 
IMPROVEMENT OF THE PERTURBATIVE 
FREE ENERGY 

The weak coupling expansion for the free energy of 
the quark-gluon plasma has been calculated through or- 
der a^^^ [il, m, , and the result yields poor conver- 
gence at the lattice QCD available temperatures T < 5Tc, 
and even at much higher temperatures. The problem is 
more involved taking into account that at order log as 
starts the contribution of non perturbative effects which 
is hard to compute [Zsj . There have also been numerous 
attempts to resum perturbation theory in order to get a 
better convergence of the result (see e.g. [1^ and refer- 
ences therein). One of the most developed techniques is 
the Hard Thermal Loop (HTL) perturbation theory, first 
proposed in [i^, . This is an efficient reorganization 
of the perturbative series which avoids the unpleasant in- 
frared singularities. We show in Fig.[T]the two- loop result 
for the trace anomaly density from HTL [28] (see also [i^ 
for a one loop computation). The HTL result becomes 



a very smooth function in the regime T > 1.13 Tj, and is 
unable to reproduce the lattice data. One of the most 
conspicuous facts is that the pressure from HTL remains 
nearly constant with a value of about 95% of that of an 
ideal gas of gluons, even for temperatures ~ 10"^ Tc- This 
means that the approach to the ideal gas is extremely 
slow, which ultimately implies that the trace anomaly 
is very small and smooth. The same feature is shared 
by the weak coupling expansion and other resummation 
techniques, and it seems to be confirmed by preliminary 
results of lattice computations at extremely high temper- 
atures J5(?|. 

Many perturbative studies based on the computation 
of the three-dimensional effective theory of QCD on the 
lattice lead once again to the same smooth behaviour for 
the trace anomaly even near Tc [Bll, [s^l . 

These problems can be understood within our frame- 
work by noting that perturbation theory, as well as HTL 
and other resummation techniques, contain only loga- 
rithms in the temperature, suggesting a mild tempera- 
ture dependence. Our discussion above shows that these 
approaches would yield a powerless contribution, Aport, 
which should ultimately be identified with the almost 
constant and vanishing oa of Eq. (j2.3p rather than with 
the full result from the lattice. Actually, the maximum 
lattice temperature T = 5Tc should be far from the 
pQCD estimate since the power correction provides the 
bulk of the full result at this temperature. For A'c = 3 
the 0{g^) corrections to E g. (|A12[) corresponds to mul- 
tiply it by (1 — 6g{n)/TT) which becomes small for 
g{^) ^ 7r/6 or ^ » IO^^Aqcd- This delayed onset of 
perturbative QCD is not a new phenomenon and takes 
place in the study of exclusive processes at high ener- 
gies where there appear collinear divergences (see e.g. 
Ref. [5^ and references therein). In order to support 
this statement, we zoom in Fig.[3]the two-loop result for 
the trace anomaly density in HTL and compare it with 
our value oa = —0.02(9) from a fit of lattice data us- 
ing Eq. (|2.3p . estimated in the continuum limit. The 
HTL result is inside the error band for temperatures 
T > 1.4(3) Tc. From a fit of = 8 lattice data we 
got OA = —0.02(4), and the HTL result is inside the 
error band for T > 1.7(4) Tc. It is appropriate to men- 
tion at this point that the estimated error of the 2-loop 
HTL start to be important below 1.5 Tc as we can see 
in Fig. [21 Moreover below 1.5 Tc the 1-loop and 2-loop 
computations in HTL start to be in disagreement [28]. 

There are in the literature other computations of the 
equation of state of gluodynamics from lattice simula- 
tions that are in qualitative agreement with Boyd et al. 
(see Ref. [lEIiSS, HH and references therein). Power 
corrections have also been obtained in the equation of 
state of full QCD [5^1 . Nonetheless we will focus in this 
paper only on pure gluodynamics for simplicity. 

The conclusion of the above discussion is that the be- 
haviour of Eq. (|2.3p clearly contradicts perturbation the- 
ory which contains no powers but only logarithms in the 
temperature. Power corrections are understood as the 
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FIG. 3: The trace anomaly density (e — SP)/T'^ as a function 
of T (in units of Tc). The 2-Ioop HTL resuh from ref. [28] is 
shown as a shaded band that corresponds to varying the scale 
/i by a factor of two around 2ttT. The horizontal lines display 
the fitted value of a a in Eq. (|2.3p . from A^r = 8 lattice data, 
a A = —0.02 ± 0.04 (error bar in green), and the continuum 
limit estimate, a a = —0.02 ± 0.09 (error bar in red). 



high energy trace of non-perturbative low energy effects. 
This kind of corrections have aheady been detected and 
analyzed in the Polyakov loop and in the heavy quark- 
antiquark free energy in terms of dimension two gluon 
condensates d, The rather successful description of 
the lattice data for the trace anomaly suggests a similar 
approach for the equation of state, that will be analyzed 
in Sec. El 

Our previous numerical analysis of lattice data sug- 
gests that the perturbative series of the free energy of a 
hot gluon plasma is rather small, not only in the high 
temperature limit, but also near the deconfined phase. 
The perturbative series has a problem of convergence for 
temperatures below 10 Tc, which could mean a lack of 
analyticity near the phase transition. 

We can improve the convergence by considering a re- 
summation of the perturbative series based on renormal- 
ization group invariance. The details of the computation 
are provided in Appendix [Cl In practice this means to 
reorganize the perturbative series and express it in terms 
of the following manifestly Renormalization Group (RG) 
invariant 



a(r) 



1 



1 



(4.1) 



The result for the pressure so obtained is a manifestly 
RG invariant series up to order 0{a^) and is given by 
Eq. (|C10|) with coefficients Eq. (|C11|) . The only unknown 
coefficient to this order is Aq, which enters at 0{a^). The 
trace anomaly involves a derivative of the pressure with 
respect to temperature, see Eq. (|2.2p . Using (|4.ip . this 
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FIG. 4: The trace anomaly density (e — 3P)^r* as a function 
of T (in units of Tc). Lattice data are from ;9] for A'^^ x Nt = 
32^ X 8. The fit using Eq. (|43)) is plotted. We also plot the 
perturbative contribution from Eq. (|4.2p , and compare it with 
the 1-loop [ii] and 2- loop '28| HTL resuhs (shaded bands). 



can be computed as 



^pcrt 



45 



-7oa 



~ 2 ^-^pcrt 



da 



(4.2) 



where we have defined TJpcrt = ^'pcrt/-Pidcai, being 
P^^c^i = {{N^ — l)7r^/45)T* the pressure of an ideal gas of 
massless gluons. This result holds modulo 0{a^^^) when 
Eq. (|G10[) is considered. 

One can try to compare the RG-invariant Apcrt, trun- 
cated to 0{dJ^'^) (included) with the available lattice 
data. In all cases we use the more accurate Nr = 8 data, 
including points in the range 1.13 < T/Tc < 4.5. To be- 
gin with we allow both Aq and Aqcd to change as free pa- 
rameters. In this form no acceptable fit is achieved. The 
best fit gives x^/dof = 1.6 with Ae = 0.33 ±3.0. It mim- 
ics the overall shape of the data but fails to reproduce 
the highest temperature data, where perturbation theory 
(PT) is expected to be more reliable. In addition, this 
best fit requires unphysical values of Aqcd ~ 338 MeV, 
which would not reproduce other gluodynamics data. To 
avoid this problem, in what follows we fix Aqcd to lie on 
the range 0.877(88) Tc, allowing for a 10% uncertainty in 
this parameter. In this case, fitting Aq in order to repro- 
duce the highest temperature data, one can verify that 
the x^/dof rapidly deteriorates as new lower temperature 
data are being included in the fit, while simultaneously 
the fitted value of Aq rapidly changes to larger negative 
values. 

Next we proceed to add the non perturbative (NP) 
term 6a(Tc/T)^ to the perturbative ones and try to re- 
produce the data using Aq and 6a as parameters. In a 
more detailed treatment one should expect some inter- 
ference between these contributions, from radiative cor- 
rections in the form of anomalous dimension, etc, in the 
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FIG. 5: Correlation ellipses corresponding to Ax^ ~ 2.3, 4.7 
(dof — 6) between the perturbative parameter Ae and the 
non perturbative parameter &a (see main text) characterizing 
the trace anomaly. The central values are obtained using the 
lattice data in the interval l.lSTc <T < 4.547^. 




1 2 3 4 5 



T/Tc 

FIG. 6: The pressure P as a function of T (in units of Tc). 
Lattice data are from ^9] for x Nr ^ 32^ x 8. We plot 
the fit of Eq. (|4.4|l , and the analogous equation with Phtl to 
2-loops [2^. We retain the same values of 6a as in Fig. 3] 
The perturbative contribution with no extra contributions is 
also displayed. 



NP term. Presently, we adopt the simplest scenario of 
additive PT and NP contributions, i.e. we consider 

A = Apert + &A y) • ^^-^^ 

This procedure yields a fairly good fit to the data, 
xVdof = 0.40, with &A = 3.18(74) and Ae = 20.0 ± 10.5 
(see Fig. 2]). This fit is only slightly worse than that 
in section [U which used a constant perturbative back- 
ground a A- The error bars in 6a and Aq are enhanced as 
they include the uncertainty in Aqcd- It is also note- 
worthy that the two parameters are highly correlated 
(see Fig. [5]), and in fact the combination Aq — 14.2 6a = 
—25.1 ± 1.3 (which has zero correlation with 6a) has a 
much smaller error. Nevertheless, we caution that these 
numbers are to be taken as indicative only. The missing 
terms in the perturbative expansion make the pertur- 
bative contribution to be little reliable at temperatures 
near the transition. For instance, if the previous fit is 
repeated using data in the range 1.43 < T/Tc < 4.5, the 
central value of 6a increases to 3.47 while Aq changes 
correspondingly. 

A remarkable fact is that the central value of Aq is 
such the Aport falls on top of the 2-loop HTL calculation. 
Actually, the value of Aq required to produce this match 
between the two perturbative calculations at T = 4.5 Tc 
is Aq — 19.0 ± 1.2 (note that this latter number does not 
depend on lattice trace anomaly data). Similar values, 
Aq = 20.0 ± 0.8, allow to reproduce the HTL result by 
Apcrt in the range 2.75 < T/Tc < 4.5 with x^/dof = 
0.67 (where we include the uncertainty from Aqcd and 
from the scale fi in HTL). While the central values of 
Aq obtained from a fit to lattice data, including the NP 



term, and to 2-loop HTL are equal (20.0), the uncertainty 
assigned from the second determination is much smaller. 
(Nevertheless, the preferred value of A^ as extracted from 
a match to HTL evolves at very high temperatures to a 
value close to zero.) 

We have also considered a fit to the lattice trace 
anomaly data using the 2-loop HTL result as pertur- 
bative background, plus the non perturbative term, in 
the range 1.24 < T/T^ < 4.5. This gives 6 = 3.62(11) 
with x^/dof = 0.72. So it also provides a good fit of the 
data, only slightly worse than that from the RG invariant 

Apcrt ■ 

In all cases the perturbative contribution is subdom- 
inant and very small above 2Tc. Also, the pressure of 
gluodynamics is very well reproduced for all tempera- 
tures above Tc, after inclusion of the non perturbative 
term. Integrating the trace anomaly, as modeled in (j4.3p . 
yields the following form 

P{T) , Ppo.,(T) bA fTcY 

where Cp is an additive constant related to the normal- 
ization of the thermodynamic quantities on the lattice. 
Note that Ppcrt(T) is known from Eq (|CTO)) . A fit of the 
lattice data for the pressure using Eq. (|4.4p is plotted in 
Fig. m We retain in the fit the same value of 6a that 
we previously obtained from the trace anomaly, and con- 
sider as free parameters Aq and Cp. This formula yields 
a good fit to the data in the regime 1 < T/Tc ^ 4.5 with 
Ae = 20± 10, Cp = -0.08 ± 0.08 and xVdof = 0.26. 
As expected, the value of Aq agrees with our previous 
determination from the trace anomaly. We also consider 
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a fit of the lattice data using the 2-loops HTL results for 
the pressure, Phtl The only free parameter in this 
case is Cp, and we get Cp = 0.006 ± 0.020, x^/dof = 0.96. 
We have increased the errors for the pressure of Ref. Q 
by a factor 2.5, and it seems to be reasonable taking 
into account the difference between lattices with differ- 
ent temporal extent N-r- 

We have also analyzed different lattice data other than 
those primarily considered above [9[, namely Ref. 
which almost reproduces and Ref. [13] which dif- 
fer from those of Ref. @ in the higher temperature re- 
gion, possibly due to systematic errors not provided in 
Ref. ^0[ . An estimate of the systematic errors based on 
a variation of lattice sizes suggests that they could be 
of the order of 0.3 in A, which explains the discrepancy 
of 0.15 between both lattice data. Thus, we have artifi- 
cially increased the errors of Ref. by a factor of ten. 
For the resummed PT-I-NP approach we get Aq — 31(7) 
and &A = 3.9(0.5) with xVdof = 0.40 for the data of 
Ref. [l| while Aq = -9(18) and &a = 1.42(1.45) with 
X^/dof = 0.12 when data from Ref. are used. Note 
that the resulting pairs (A6,6a) fall reasonably close to 
the line Ag- 14.2 6a = -25.1 ±1.3. The HTL-^NP yields 
&A = 3.44(0.09) with xVdof = 2.6 for Ref. [l|, and 
6a = 4.05(0.15) with xVdof = 0.44 for Ref. With 
the provisos spelled out previously, it is quite reasonable 
to conclude that other lattice data sets confirm our re- 
sults, with larger uncertainties. 

To finish this section we note that the trace anomaly 
as obtained from Aport or from HTL does not have a 
definite sign, but this quantity becomes positive definite 
for all temperatures above upon inclusion of the non 
perturbative term bA{Tc/T)'^, for values of 6a fitting the 
data. Curiously, the fitted value Aq = 20 is roughly the 
threshold value above which A would fail to be positive 
definite (the region with negative values being located 
around T = 50 T^). 



V. DIMENSION TWO CONDENSATE AND 
TRACE ANOMALY 

We can make a first approximation to the computation 
of the trace anomaly using the non perturbative model in- 
troduced in There, the inverse temperature power 
corrections observed in lattice data for the Polyakov loop 
[63 | and the gg-potential [l^ , are accounted for by means 
of a non perturbative contribution to the zeroth compo- 
nent gluon propagator of the dimensionally reduced the- 
ory. Namely, 



^00 (k) 



(k) 



rjNP 



(k), 



(5.1) 



non perturbative term in the gluon propagator induces 
a corresponding contribution to the dimension two gluon 
condensate: 



,NP 



l)Tm% 



(5.3) 



Being the trace anomaly in giuodynamics given by 
Eq. (|2.ip , the problem can be addressed through the com- 
putation of the vacuum expectation value of the squared 
field strength tensor ((GJ^^)^). 

In perturbation theory (G^^^) starts at 0{g^), cf. 
(IA12p . This comes from integration of "hard modes", 
i.e., modes with momentum scale 2ttT [6^. The con- 
tribution of soft modes, described by the dimensionally 
reduced Lagrangian, starts at 0{g^). These perturbative 
contributions break scale invariance through radiative 
corrections and so depend logarithmically on the temper- 
ature, after extracting the canonical factor T*. Power 
corrections require a stronger breaking of scale invari- 
ance. We will assume that this does not happen for the 
hard modes since condensates are low energy phenom- 
ena. Therefore, we assume that the non perturbative 
scale breaking enters through the parameter uiq in (|5.2p 
corresponding to the temporal gluon propagator. Just 
by dimensional counting, this parameter will produce the 
required power correction, ~ in the observable A. 

The ultraviolet divergent momentum integrals in three 
dimensions are dealt with by dimensional regularization. 
On the other hand no non perturbative scale breaking is 
assumed for spatial gluons. As these gluons fail to have a 
Debye mass, a non perturbative term of the type m^^/fc^ 
in the propagator of spatial gluons would yield a vanish- 
ing contribution within dimensional regularization. This 
approach yields 



2nNP 
00 



(5.4) 



= -T(7V2 - 1) Amz^m^j. 



The negative sign is a consequence of the renormaliza- 
tion.^ Making use of (|5.3p . this result can be expressed 
in a more conveniently form in terms of the dimension 
two gluon condensate, and finally gives 



(e - 3P)np = -3 
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)(^0,a 



NP 



(5.5) 



where 



The behaviour is T^, in contrast to the perturbative 
behaviour ~ T'*. It is also noteworthy that (e — 3P)np 
satisfies the correct large Nc counting, being 0{N^). 



<(k) = 



k2 



i?oT(k) = 



(k^ -I- m|,)2 



(5.2) 



Here mu is the Debye mass, which scales like T, and 
ma is a temperature independent mass parameter. The 



As advertised, an 0{g'^) contribution, namely, proportional to 
m|j, is derived from the analogous computation using the per- 
turbative component of the gluon propagator. 
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The coefficient &a in the parameterization (|2.3p is to 
be identified with the estimate of Eq. (|5.5p ^ 

bAT!^-S^'^{AlJ^^. (5.6) 
9 

Therefore, since I3{g) < 0, we have a positive &a- If 
we consider the perturbative value of the beta function 
Pig) ~ 5^ + 0{g^), the r.h.s. of Eq. ([511) shows a fac- 
tor in addition to the dimension two gluon condensate 
(7^(Aq (j)^^. So the fit of the trace anomaly data is sensi- 
tive to the value of g. Fortunately, this quantity has only 
a smooth dependence on T. For the Polyakov loop the 
sensitivity in g is only through the perturbative terms, 
which are much smaller than the NP ones. When we 
consider the perturbative value gp to 2-loop, gp ^ 1.41, 
we get from the fit of the trace anomaly 

g'iKaf'' - (0.72(3) GeV)^ (5.7) 

Here we have used the continuum limit estimate 6a = 
3.28(27) of section ini based on a constant perturbative 
background. For the present purposes we consider this 
estimate sufficiently accurate and less subject to details 
related to the precise determination of the perturbative 
background. The precise number for 5a varies by using 
any of the determinations of 6a discussed in section IIVI 
but they are still compatible with each other. 

The value quoted in Eq. (|5.7p is a factor 1.5 smaller 
than the value obtained from a fit of the Polyakov loop 
(0.84(6) GeV)^ [1] and heavy qq singlet free energy 
(0.90(5) GeV)^ This disagreement could be explained 
in part on the basis of certain ambiguity of g in the non 
perturbative regime. We show in Ref. p that a better fit 
of the Polyakov loop and heavy quark free energy lattice 
data in the regime Tc < T < ATc is obtained for a value 
of g slightly smaller than gp, i.e. g — 1.26 — 1.41. Taking 
this, we get 

g'iK,)''^ = (0.77(6) GeV)2 

= (2.84 ± 0.21 Tc)2 , (5.8) 

in better agreement with determinations of the conden- 
sate from previous observables, see Table H] This value 
follows from the continuum limit estimate of 6a, and its 
discrepancy with the corresponding value from a fit of 
Nr = 8 lattice data is 5%, which is much smaller than 
the error. The same happens with the lattice data of the 
Polyakov loop [3] . 



^ It may look unnatural that the rather complicated r.h.s. of 
this formula yields b^, when Pisarski's argument, section IlIII 
assigned to this quantity an almost geometrical value, (A'^^ — 
l)27r^/45. However, this is not so. In the present model, the ge- 
ometrical value of 6a together with the r.h.s. provide rather an 
estimate of Tc, identified as the temperature where the pressure 
becomes negligible. 



Observable /(^o,a)^^ 

Polyakov loop [4] (3.11 ± 0.22 

Heavy qq free energy [6] (3.33 ± 0.19 Tcf 

Trace Anomaly (2.84 ± 0.21 T^)^ 



TABLE I: Values of the dimension two gluon condensate from 
a fit of several observables in the deconfined phase of gluody- 
namics: Polyakov loop, singlet free energy of heavy quark- 
antiquark and trace anomaly. Values are in units of Tc. We 
show the fit for lattice data with Nt = 8 for the heavy qq free 
energy, and the continuum limit estimate for the others. Error 
in last line takes into account an indeterminate value of the 
coupling constant g — 1.26— 1.41, being the highest value the 
perturbative gp up to 2-loop at T = 2Tc. The critical tem- 
perature in gluodynamics is taken as Tc — 270 ± 2 MeV [6^ . 



In a recent paper the electric-magnetic asymmetry of 
the condensate defined as 

AMT)=g^{AlJ-^g^{Al,), (5.9) 

is computed in SU(2) Yang-Mills theory [63|- Obviously 
this quantity vanishes at zero temperature by Euclidean 
symmetry. The dimension two condensate has dimension 
of mass squared and so it would vanish in a perturba- 
tive calculation at zero temperature. At finite temper- 
ature instead (^oa)^ scales as modulo slowly vary- 
ing radiative corrections. On the other hand our model 
assumes that the non perturbative part of the conden- 
sate {Aq a)^^ is temperature independent (modulo ra- 
diative corrections) at least in the regime slightly above 
the phase transition and beyond This is exactly the be- 
haviour that is obtained in [67J for the asymmetry in 
this regime. If we naively scale our determination for 
g^(AQjj)^^ with the number of gluons 3/8, we get in 
SU(2)'the estimate g^iAl J^^ = {l.9{2)Ta)^. This value 
can be compared with the temperature independent con- 
tribution to A^°i = (0.894(14) Tc)2 of [67], leaving 
some room for a magnetic contribution in our model. In 
any case, the very existence of the dimension two con- 
densate seems to be confirmed by the work of ref. [67| . 



VI. CONCLUSIONS 

For quite a long time it has been believed that the high- 
est temperature measured in lattice calculations Tmax ^ 
5 Tc could be matched smoothly to perturbative calcu- 
lations. Assigning the standard MS conversion factor to 
a momentum scale fi ^ 2'kT and taking Tc ~ 0.6y^ 
yields /Xmax ~ 1.5 GeV; a high scale. Actually, there 
have been many perturbative studies that unsuccessfully 
tried to reproduce the lattice data for the trace anomaly. 
As we have shown in this work, the trace anomaly den- 
sity shows unmistakable traces of power corrections in 
the inverse temperature in the region slightly above the 
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phase transition which may be accommodated by dimen- 
sion two gluon condensate. This finding parallels a simi- 
lar analysis for the Polyakov loop and the values for the 
dimension two condensate is in satisfactory agreement 
with the value found here. The 20% discrepancy may 
be due to many different sources such as anomalous di- 
mension effects. It is important to emphasize that the 
non-triviality of such a numerical consistency can only 
be addressed in a calculation where all quantities can be 
estimated simultaneously. In addition, it would be wor- 
risome if the dimension two condensate associated with 
the power correction would differ by an order of mag- 
nitude. The fact that this is by far not the case sug- 
gests searching for other features where the condensate 
might generate further genuinely non-perturbative ther- 
mal power corrections. This of course applies not only to 
gluodynamics but also to the case of full QCD with ac- 
tive dynamical quarks and in the presence of a chemical 
potential. A possible guideline for future studies would 
be a simple minded extension of the power corrections 
ideas discussed in the present work. Actually, most re- 
cent QCD calculations already confirm the ubiquitous 
1 /T^ behaviour of the trace anomaly density [sl] . 

On a more methodological level, let us remind that 
at high energies and zero temperature momentum power 
corrections are widely accepted and encode via conden- 
sates non-perturbative effects. For this situation the full 
machinery of QCD sum rules has been extensively de- 
veloped and applied with recurrent success. While there 
is no reason why condensates should not be present at 
high temperatures as genuinely non-perturbative effects 
even after deconfinement, temperature power corrections 
above the phase transition are much less common and a 
more complete technology than that used in the present 
work would be most useful. In addition to the standard 
finite temperature complications one must also face the 
well known ambiguity on clear separation between per- 
turbative effects from existing non-perturbative features. 
There is no doubt that this intricate theoretical problem 
should necessarily be addressed at some stage in the near 
future. 

In spite of the overall numerical consistency found for 
the dimension two condensate, the physics underlying the 
generation of such a condensate, and ultimately the very 
appearance of thermal power corrections remains uncer- 
tain. Actually, the elusiveness of a theoretical determina- 
tion of the dimension two condensate from first principles 
refiects our current ignorance and remains a challenging 
bottleneck to the understanding of strongly interacting 
quark gluon plasma in non Abelian gauge theories. 
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APPENDIX A: TRACE ANOMALY AT FINITE 
TEMPERATURE 

In this appendix we review the derivation of the trace 
anomaly at finite temperature [H, [12, IH, [2^ [2^. In 
the Euclidean path integral representation, the partition 
function for gluodynamics reads 



Z = / VA^^a exp 



(Al) 



where the field strength tensor is given by (5^,^ = d^A^ — 
d^Afj, — i[Af_i,,A^] and (G^^,)^ > in Euchdean space. 
Gauge fixing terms will be considered in the giuonic mea- 

The 
This 



sure DAi^^^a, with periodic boundary conditions 



canonical gluon fields are given by 
provides the relation 



dlogZ _ 1 
dg " V 



j_u 



{{G%f) (A2) 



where in the last equality translational invariance has 
been used. 

From the standard thermodynamic relations, the free 
energy, the pressure and energy density are given by 



F = ~PV = -TlogZ 

_ E dlogZ 
^ ^ V ^ Y dT 



as well as the relation 



3P = T' 



d (P 



dT 



(A3) 
(A4) 

(A5) 



The dimensionless quantity P/T^ is a function of 
T/Aqcd (or rather Agiuodynamics in our case) and in prin- 
ciple it can be expressed to all orders in terms of the 
coupling constant and the scale: 



^ = f{g{^^)^og{^l/27^T)). 



(A6) 



The right-hand side is actually independent of the scale 
fi. This allows to compensate a variation in T with a 
variation in fi and then in g (i.e., in Aqcd) 



d 



P 



dg d f P 



dlogT J dlogfidg\T^ 



(AT) 
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Using now 



P 

y4 



\ogZ_ 



(A8) 



with (|A2p and (jA5P yields finally the trace anomaly, 

Pig) 



3P 



25 



where we have introduced the beta function 

dg 



Pig) = 1^ 



djjL 



(A9) 



(AlO) 



In this work we focus on the evaluation of the trace 
anomaly density above the phase transition, which in 
practice means momenta above ~ 1 GeV if the cus- 
tomary MS conversion factor fi — 2ttT is adopted. In 
perturbation theory the beta function reads 



Pig)^-P^g^ + 0{g'') 



(All) 



where Pa = lliVc/(487r2). 

The result for the trace anomaly density to two loops 
is [si 



e-3P 



72 



/?o/(T) + 0(g'^) (A12) 



where l/g'^{fi) — /Jq log(/i^/AQQj-,) at leading order. Re- 
markably this perturbative result implies a negative ex- 
pectation value of the squared field strength tensor; this 
quantity being positive prior to renormalization. Note 
the ambiguity in the (truncated) perturbative result, 
since generally one has both the temperature T and the 
MS-renormalization scale /z, for which it is usually taken 
in the literature the reasonable but arbitrary choice 
/i ^ 27rT. This issue is analyzed in section IIVI Higher 
order corrections including up to g^ log g can be traced 
from [51j . Unfortunately severe infrared problems in the 
perturbative expansion yield poor convergence at tem- 
peratures available to lattice QCD calculations, T < 5Tc. 



APPENDIX B: THE GLUEBALL HAGEDORN 
SPECTRUM 

The known glueball spectrum as obtained from the lat- 
tice ends at a maximum mass value Mmax- Although this 
is somewhat tangential to the subject of the present pa- 
per it is interesting to analyze the effects due to a possible 
contribution of higher mass states. The density of states 
is defined as 



p(Af)=^g,5(M-M,) 



dN{M) 
dM ' 



(Bl) 



where N{M) is the cumulative number of glueball states 
NiM) = J2 g^^iM - Mi). (B2) 




Mm 



FIG. 7: Logarithmic plot of the cumulative number of glue- 
ball states obtained from the lattice as a function of the 
mass in units of the Sommer parameter ro compared to the 
exponential spectrum fit N(M) = Ae'^^'^" (straight line) 
with Th = 2.1/ro. The smoothed cumulative number with 
vqAAI — 0.5 is also shown. 



Here we conventionally take 0(0) = 1/2. A logarithmic 
plot of N{M) for the glueball spectrum is presented in 
Fig. [7] where a rough straight line can be envisaged with 
the exception of low lying states and the highest states, 
presumably due to boundary effects.'^ We fit the function 



N{M) = Ae"/^" 



(B3) 



by minimizing 



J2 (logAf(M,) - log(^e^^'/^«))' . (B4) 



This gives roTn = 2.1 or equivalently Th = 2.8 T^. A 
weak point of the present determination is the inabil- 
ity to provide a reliable error estimate for the former 
fit, as there are fluctuations in the spectrum that our 
smooth fitting function can never account for. To con- 
sider this possibility we also represent the smoothed cu- 
mulative number 



{NutiM)) =Y,g^ f-tan-1 



M - M, 
AM 



l) (B5) 



which mimics a Breit-Wigner finite width in the density 
of states. Fluctuations in the spectrum are washed out 
already with a common smoothing tqAM = 0.5, a factor 
only 5 times larger than the uncertainty in the lowest 
glueball mass. As we can see, the smoothing tends to 
confirm the fit, suggesting the robustness of the analysis. 



^ For the skeptical reader we note that recent updates of the par- 
ticle spectrum provide a similar pattern for both mesons and 



baryons separately 1681 16 
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APPENDIX C: PERTURBATIVE EXPANSION 
OF THE FREE ENERGY 

We compute in this appendix the perturbative ex- 
pansion of the free energy of a hot gluon plasma using 
a renormahzation group (RG) improvement of the se- 
ries. The pressure of gluodynamics has been calculated 
within the weak coupling expansion through 0(g^ logg) 
in Ref . 0, |4^, [H, |53, HE Hi) US HJ , up to an unknown 
constant. The result can be written in the form 



pert 



1 + + a^a' 

„2 



3/2 



04 



.5/2 



&4log(2:^) +C4loga 
a5 + 65log(^) 

^elog^^:^) +celoga 
+d6 log {■^) log" + 66 log^ ( 



2ttT 



(CI) 



where we have defined iiTport = -PpertZ-Pideai, being 
^'idoai = ((-^c ~ l)7r^/45)r'' the pressure of an ideal gas of 
massless gluons. In this formula a = = (7^(/z)/47r is 
the running coupling constant. All coefficients appearing 
in (jCip are known except ag , which cannot be computed 
within the perturbative scheme owing to infrared diver- 
gencies [6^ . 

The trace anomaly follows straightforwardly from the 
expression of {{N^ - l)7^V45)i^pert by applying Td/dT 
(the derivative affects only the explicit T dependence). 

We can see in Eq. (jCip that Ppcrt{T) is computed in 
the weak coupling limit as an expansion in powers of a(/i) 
with coefficients which depend on log(/i/27rT). However 
the dependence in ^ disappears when all orders are added 



d 

/^^^^port — 0. 



(C2) 



This is because not all coefficients a^, &i, Ci, are indepen- 
dent. Indeed 



= 7(a) 



with 



7(0;) = —a^ (70 -I- 71a + 72^^ + • • •) 



(C3) 



(C4) 



and 70 = f , 71 = ^ and 72 = ^ for - 3. 
When this is applied to Eq. (jCip it implies 

hi = 7oa2, 

, 3 

h = 7o(2a4 + C4) + 7ia2, 

da = 270C4, 

66 = 70^2 • (C5) 



Of course, RG invariance is spoiled if only a finite number 
of terms is retained. For the truncated series to be useful, 
a{^) must be small, i.e., /i 3> Aqcd, with fixed 2'kT / 
which in turn requires T ^ Aqcd- 

We can go further and try to reorganize the series of 
Hpcrt in such a way that it appears explicitly RG in- 
variant. To this end we need to identify some invariant 
quantity and make the expansion around it. The solution 
of Eq. (jCSp can be written as 



da 



(mo) 

7(a) 



log 



(C6) 



Defining, as usual, Aqcd by the condition a(AQCD) = 00, 
implies 



(m) 



da 
7(a) 



log 



A 



2ttT \ 



QCD 



(C7) 



This allows to express a(/i), appearing in the left-hand 
side, in terms of log(AQCD/27rr) and log(^/27rT). For 
convenience, we introduce the quantity 

a(T) = l T^TTT' (C8) 



^0 1O£ 



( Jji. 



which is manifestly RG invariant and behaves as a(/i) 
when /i = 2ttT 00. From jCTI and ([Ci]) then 



;r2 



Q;(/i) = a + a 



log5-7ofog(^) 



+a^ 



72 
70 



^'+4loga- 



-^^l°s(^)+7o^fog^( 
-271 log ci log 



4 log' « 

7o 



2tiTI 



2ttT 



(C9) 



By introducing this perturbative series of a(/i), in 
Eq. (jCip one obtains 

A2a + Asa^/^ 
[A4 + B4 log a] 



pert 



-52 



-5^/2 + log 5] 
-a' 



[^6 + ^6 loga + Celog^ a] 
+0(a^/2)_ (c^o) 

All coefficients here are known except Ag. In fact'' 
A2 = a2, 



^ The values for the coefficients are A2 = -15/(47r), A3 = 30/7r3/2, 
A4, = 135(3.51 - log7r)/(27r2), B4, = 5175/(887r2), A5 = -3.23 * 
495/(27rS/2), B5 = 2295/(227r5/2), fig = 2295/(19367r3)(115 -|- 
264(3.51 - logTr)) - 1134.8/7r3, Cs = 566865/(19367r3). 
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^3 = as, 

A4 = 04, 

B4 = — 02 + C4, 

70 

A5 = 05, 

270 

72 li\ 

2 + ae, 

70 7o/ 

7i 7i 7i 
Be = -,02 + 2 — a4 H Ci + c%, 

7o 70 70 

Oe = — ^a2 + 2 — C4. 
7o 



70 



(Cll) 



Hpcrt is now written in a manifestly RG invariant way 
as a function of T/Aqcd- In principle, the parameter 
Aqcd itself can be extracted from the pressure at suffi- 



ciently high temperature: 

1 



log A, 



QCD 



lim 



pert 



1 



i2 



70^2 T^oo 
--1J2 --B4l0gt- ^4^, 



A2 
t 



(C12) 



where t = 1/(70 log(27rT)) (and Aqcd results in the same 
units used for T in t). In practice, the pressure is not 
known to the required precision to such high tempera- 
tures and the parameter Aqcd is obtained from other 
determinations. 

The trace anomaly follows immediately from the RG 
expression noting that da/dlogT = —joa^, 



*pert 



45 



~2 "-f^pert 

7oQ —pr- 
da 



(C13) 



known modulo 0{a^^'^). 
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